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ABSTRACT 


Expendable Bathythermograph (XBT) temperature depth profile data 
from 17 locations in the north Pacific Ocean has been analyzed to extract 
information concerning the thermal microstructure. The microstructure 
information is then examined to determine if the distribution can be 
associated with the climatic regions of the north Pacific Ocean proposed 
by Tully [ 1964 ]. 

A modified binomial filter used by Roden [ 1971 ] to analyze STD data 
was applied to find a mean temperature series, A temperature per- 

ixbation series results when the mean is subtracted from the original 
XBT series. Two resulting parameters, average mean ternperature 
perturbation vertical thickness scale and average maximum amplitude 
of the temperature perturbations, can be associated with the micro- 
SEPuectuUre. 

The average mean temperature perturbation vertical thickness scale 
was found to increase with latitude from 20°N to 55°N, and as the thick- 
ness scale increased, the maximum amplitude of the perturbations was 


found to decrease. 
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I, INTRODUCTION 


propagation paths for acousticel energy. The gross temperature structure 
is used for ray path calculations to predict transmission characteristics 
[Cohen and Weinberg, 1969]. However, sound measurements at one location 


reveal a fluctuating accustic pressure level {Sheehy, 1950] attributed to 


variations in the fine temperature and salinity structure that cause small 


but important changes in the propagation paths. At higher frequencies | 
‘the microstructure becomes increasingly significant. A theory for acoustic 


wave transmission has been developed for a medium containing thermal 





inhomogeneities of a thermal nature that are Spherical or spheroidal 


"patches" randomly distributed in space [For example, Liebermann, 1951; 


Mintzer, 1953a, 1953b, 1954]. The measured sound fluctuations have been 





patch" distribution assumption. A layered microstructure with layers 


of small gradient a few to a few tens of meters thick separated by sheets 


of high gradient a few to a few tens of centimeters thick ‘nd with 
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horizontal dimensions up to thousands of meters have been found in many 
regions of the oceans [Cooper, 1967; Cooper and stone 1968; Neshyba, 
Neal, and Denner , 1969; Osborn, 1969; Roden, 1971; Simpson, 1971; 
Stommel and Fedorov, 1967; and Woods, 1968]. Several mechanisms have 
been proposed for generating the layered microstructure, and it may be 
that a different mechanism is dominant for different ocean regions. 
This layered microstructure concept requires a re-examination of the 
present theory for the transmission of acoustic energy in a randomly 
inhomogeneous medium, 


In 1957, Whitmarsh, Skudrzyk, and Urick concluded a paper on the 


scattering of sound in the sea with the following comment: ".... the 





Some 11 years later, 





during a meeting on Turbulence In The Ocean held at The University of 
British Columbia, R. W. Stewart "pointed out that a case existed for 
another degree of exploration of the ocean. While in the past, widespread 
surveys have sought to measure the temperature and salinity fields, we 


must now concentrate on a widespread search for fluctuations in the 


1 Whitmarsh, D.C. , Skudrzyk, E. , and Urick, R.J., "Forward Scattering 
of Sound in the Sea and Its Correlation with Temperature Microstruc- 
tume, |. cous. Soc. Am., Vv. 29, p. 1141, October 1957. 
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temperature and salinity fields. * The only data that exists with such 
broad coverage is bathythermograph data. Magruder [1970] suggests that 
Expendable Bathythermograph (XBT) data might be useful. 

The purpose of this thesis is to see if information about the distri- 
bution of microstructure can be determined from XBT traces. If the 
microstructure distribution can be determined from an XBT trace, then 
using some of this available data, can the distribution be related to the 
climatic regions of the ocean for use as a rational basis for assessing 
the acoustic variability in these regions? Finally, what operational 
techniques can be developed using knowledge of the microstructure deter- 
mined from the XBT trace that would be of benefit to the Navy in the 
operating environment? 

Tens of thousand of XBT traces are on file at the National Oceano- 
graphic Data Center (NODC), with the number increasing each day. Under 
the Fleet Numerical Weather Central (FNWC), Cooperative Oceanographic 
Observation Program, XBT traces are also available from commercial 
carrier trans-atlantic and trans-pacific routes. While the distribution 
is by no means uniform in space or time, it is the only source of micro- 


structure data covering the world oceans. 


2 Osborn, T.R. , Report of the Symposium on Turbulence in the Ocean, 
Institute of Oceanography, University of British Columbia, Vancouver, 


British Columbia, p. 20, July 1970. 





lie THEORY OF ACOUSTIC PROPAGATION IN A 
RANDOMLY INHOMOGENEOUS MEDIUM 


‘Sheehy [1950] published data showing that the magnitude of the 


acoustic pressure maximum amplitude fluctuations (Coefficient of 
Variation) for 24 kHz direct path transmissions increased proportional 
to the square root of the range. However, it should be noted from 


Figure 1 that the scatter of data is rather large. 
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Figure 1. Coefficients of Variation of the amplitudes of direct 
signals plotted against the range to the receiver. Sound frequency, 
24 kHz [after Sheehy, 1950}. 
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Bergmann [1946] using the application of ray optics to investigate 
the propagation of radiation in a medium with random inhomogeneities , 
derived an approximate formula which showed that the signal level fluctua- 
tions increased in proportion to the threezhalves power of the range. He 
pointed out the variance with Sheehy's results, but offered no explanation 
for the differences. 

Liebermann [1951], using an autocorrelation function calculated from 
measured horizontal small-scale temperature inhomogeneities , showed 
that the refractive effects of the horizontal inhomogeneities produce a 
local warping (focusing or defocusing of sound) that causes the acoustic 
intensity to fluctuate. He found that refraction should only occur if the 
frequency is high enough to make the dimensions of the inhomogeneities 
greatly exceed the wave length of the radiation. 

Mintzer [1953a] derived a first order approximation for the acoustic 
pressure amplitude variation of an acoustic pulse traveling in a medium 
with small random variations in the refractive index that are a result of 
random temperature variations. Assuming that the pressure amplitude 
fluctuations are small enough to allow the use of a single scattering 
approximation and that the refractive index at a point only changes during 
the time between pulses and not while the pulse is passing the point, he 
solved the Born approximation to the wave equation for the pressure 
amplitude of the acoustic pulse. Then using the Coefficient of Variation, 


V, defined as: 
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BPD = 11 
KEI ° 


2 , where 


V 


<« >= ensemble average 


P =mean presure amplitude 


and by averaging the pressure amplitude over the pulse length instead of 
using the maximum value of the pressure amplitude as Sheehy did with his 
experiment:al data, Muintzer's theoretical analysis yields the following 


for ranges large compared with the correlation distance: 


2 
ves 2ko ay | N(o)do , where 
xe) 


N( 9)= autocorrelation function of the refraction 
index variation taken along the line between 


source and receiver 


Kao wave number of the sound 
r= range from source to receiver 
Q= rms refractive index variation 


Using a Gaussian correlation function the coefficient of variation 


becomes: 


x ,2 2 & 3} 
Va ( wWoekG eo alee siwhere a =ieorrelation 


HS 





This theoretical result is in agreement with Sheehy's experimental 
results [Figure 1]. 

Later Mintzer [1953b] determined over what range the first assump- 
tion (use of a single scatterer approximation) was valid. He found that 
for large ranges compared with the acoustic wave length and the inhomo- 
geneity size, and a short pulse length compared to the time required for 


a refractive index change: 


(1) The Coefficient of Variation of the pressure 


amplitude is: 


a ee 
*k, @ ar, ka>>Il 


y2 


7 
x 3 Z 


2 5a PAN 
Vo = 1m°kg a air,ka<<l 


(2) The single-scattering assumption used 1s valid 


OX: 
2. 3 
Kea een, |) 


Mintzer also made a comparison with ray theory and found that fora 


frequency high enough to satisfy: . 
2 
r/k a << 1 


the intensity coefficient of variation computed by wave theory and ray 


theory are the same. 
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In a third paper, Mintzer [1954] considered the second basic assump- 
tion, that the refraction index varies only during the time between pulses. 
He found that "if the pulse length is much less than the time constant of 
the refractive index correlation function, the coefficient of variation is 
independent of the pulse length, and the pulse-correlation function is 
found to be identical with the time-dependent part of the refractive index 


correlation function. "3 


| Skudrzyk [ 1957] discussed the pressure amplitude fluctuations caused 
by water inhomogeneities considered to be spherical "patches" that are 
a | aco’ > wat He showed that 
in the vicinity of the receiver pressure amplitude fluctuations are caused 
by focusing or defocusing from the lens action of the patches, and beyond 
a certain range the fluctuations are caused only by an interference field 
which results from deflected rays arriving at a point with phase shifts 
due to unequal paths of travel. He also pointed out that the interference 
component of scattering is independent of the shape of the patches, while 
the focusing component can be significantly affected by the shape. For 
example, a spheroidal patch having a vertical axis less than the horizontal 


axis will have a larger focusing component than a spherical patch. At low 


frequencies (below 20 kHz} interference scattering dominates, while at 


3 Mintzer, D. , "Wave Propagation in a Randomly Inhomogeneous 
Mecimn sell], JawAcousaiSoc. Am., v.26, p. 186, March 1954. 
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higher frequencies the interference component may be so small in 
comparison with the focusing component that at ranges greater than the 
focusing range the Coefficient of Variation is almost constant. The 
focusing range ( Yr, ) is the range beyond which the focusing does not 
change and additional pressure amplitude fluctuations are caused by 
interference only. Skudrzyk found the focusing range to be a function 
of frequency and path size. 

k = acoustic wave number 
r =ka? , where 

a = effective radius of patch 
This result is pictured graphically in Figure 2. If the sonar frequency 
and the patch size are known, then the focusing range can be determined. 
(For example, a 24 kHz sonar has a focusing range of 200 meters, if the 


# ae SER 
LAMY 
Ky Ne 


OAD 


: LC Ze 
AKAM LV 


an $000 





vw 


Figure 2. Focusing Range sal as a Functic 1 of Frequency 
(kHz) and Patch Size (meters) [Skudrzyk, 1957]. 
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Skudrzyk shows that the Coefficient of Variation of the pressure 


amplitude due to the focusing effect, V¢, is: 


2 3 
Va = @ (= where a is the relative root-mean- 
a 


Square deviation of sound velocity 
while the Coefficient of Variation due to the interference effect is: 


% 


us 
2. = 5 k 2% (eae) 


Vi 


Beyond the focusing range, scattering is due to both focusing and inter- 


ference so the Coefficient of Variation becomes: 


2 2. 2 
Vp Veh Von 


and for spherical patches: 


ve = a? kK? =) Ge eee 


which is Mintzer's [1953b] result except for a factor of 1/2. If the 
patches are spheroidal rather than spherical, then the focusing may be 
increased significantly. Whitmarsh, Skudrzyk, and Urick [1957] found 


the focusing to increase by a factor 


Vp Gpheroidal) _ ee ee 


Vp (Spherical) 16 


X is the ratio of vertical to horizontal axis of the 
patch 


Lo 





Figure 3 depicts this ratio graphically. The interference component of 


scattering is virtually independent of the vertical axis dimension. 


0 


V,/V, (spherical) 


r 


0.1 LO 10 


Figure 3. Increase in focusing due to spheroidal shape of a 
patch as a function of the ratio of the lengths of the vertical 


(b) to horizontal (a) axis [Whitmarsh, Skudrzyk, and Urick, 
1957]. 


Since published experimental results pertained only to one frequency 
and did not include temperature and acoustical fluctuation measurements 
made as one observation, Whitmarsh, Skudrzyk, and Urick [1957] conducted 
a series of experiments using 25. kHz and 60 kHz transmissions to provide 


additional data to test the existing theories for scattering of acoustical 
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energy in the sea. They used a gradiometer bar suspended from a buoy 

(a 12 foot long bar with 9 thermistor beads) to record temperature 
gradients as small as 0.001°C and to record the temperature fluctuations 
as a function of time. They point out that all the existing theories for 
scattering of sound in an inhomogeneous medium are based on the assump- 


tion of a spherical patch thermal inhomogmeity (one such patch is illustra- 





ted in Figure 4). — 





mm > 
oe 
i 


I 
| 


memes m~., = _ 
Les o/a res F 
C koe (2 
Ve2Vi/2 a (r/0)°/2 Vy (r4AV2 Jka Var=r,) 
FOCUSSING RANGE INTERFERENCE RANGE 


Figure 4, Illustration of the various types of acoustic scattering 
due to a thermal patch in the sea [Whitmarsh, Skudrzyk, and Urick, 
S37 |: 
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_ geneous turbulence, see Figure 5. However, even replacing the correlation 
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Figure 5. Horizontal temperature distributions as measured in the 
sea [after Whitmarsh, Skudrzyk, and Urick, 1957]. 


function by a new function representative of the Kolmogorov distribution, 
the theories for sound transmission can still be considered valid. 
Whitmarsh, Skudrzyk, and Urick concluded from their experiments that 
the present theories are capable of predicting the approximate variation 
in the pressure amplitude to be expected for acoustical transmissions in 


an inhomogeneous medium. 
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‘Sagar [1960] analyzing the experimental data he collected off New 
Zealand suggests that sound intensity fluctuations which obey the Mintzer _ 





night. He also found and described the water to have a "sandwich" or 


layered structure rather than the classical "patch" structure. 


Zo 





line RECENT MEASUREMENTS OF MICROSTRUCTURE 


A. REVIEW 

Stommel and Fedorov [1967] used a salinity, temperature, and depth 
recorder (STD) to take measurements of the fine structure (micro- 
structure) in the main thermocline off Timor in the Arafura Sea and 
Mindanao in the Philippine Sea. They were investigating reports of fine 
structure which began to accumulate with the advent of an instrument 
like the STD, which continucusly records temperature and salinity as a 


function of depth. They found the vertical microstructure in the main 





thermocline to consist o: 


(2 to 5 meters thick) of high temperature gradient. As an example, at 


one station a mean thermal gradient of 0.040°C/meter was found, while 





in a homogeneous layer the gradient was greater than 0.35°C/meter. By 
tracing some of the well defined laminae in a horizontal direction, they 
found the laminae extended from 2 to 20 kilometers. Though the STD 
resolution was only 2 meters, Stommel and Fedorov felt the micro- 
structure existed on a smaller scale and that further investigation would 
be of interest. 

Cooper and Stommel [1968] made microstructure measurements in the 
main thermocline off Bermuda using an STD. They too found a layered 


structure with the homogeneous layers from 3 to 5 meters thick and the 
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high gradient layers between 10 and 15 meters thick. In the layers of 
high gradient, temperature changes as high as 0. 50°C were found. The 
horizontal extent of the layers was between 400 and 1000 meters. This 
type structure was present out to a distance of 120 kilometers from 
Bermuda. Based on their measurements, Cooper and Stommel suggested 
that most of the Sargasso Sea exhibited the same characteristics in the 
main thermocline. 

Woods [1968a] conducted an extensive study of the microstructure 
in the summer thermocline of the Mediterranean Sea. A temperature 
eradient meter was used to obtain measurements of temperature gradient 
while being lowered and temperature while being recovered. The instru- 
ment consisted of two thermistors mounted 50 centimeters apart ona 
vertical fiber glass rod. Only one of the thermistors was used for the 
temperature measurements. The instrument had an 0.01°C accuracy and 
a depth limitation of 50 meters. Woods found layers with a moderate 
temperature gradient that were thick in comparison with alternate layers 
(sheets) of high temperature gradient. Figure 6 is a plot of his results 
for one sounding. High gradient sheets can be seen at depths of 10, 19, 
24, 32, 41, and 50 meters. By injecting dye from suspended packets into 
the different layers, Woods was able to associate laminar flow to thin 
high gradient layers and a weak turbulent flow to the thick homogeneous 
layers. The results of these two different type flows is a strong sheer 


across the high gradient sheets. 
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Figure 6. Temperature and temperature-gradient sounding 
for 0630 24 August 1966 in the Mediterranean Sea [Woods, 1968a]. 


The next year Woods [1968b] conducted another investigation of 
vertical heat flow through the upper ocean in the Mediterranean Sea off 
Malta. Again the sheets in the thermocline were found as can be seen 


in Figure 7, 
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Figure 7. Temperature and temperature-gradient soundings 


for Maltese waters, 14 September 1967 [Woods, 1968b]. 


Lovett [1968] took temperature profiles in the California current 


off Baja California using a sound, velocity, temperature, and pressure 
recorder (SVTP). He too found a definite layered structure in all the 
temperature profiles. He suggests that the sharp gradients found in 


the step-like structure may have an effect on acoustic scattering, 


especially so at high frequencies. 
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Osborn [1969] using a free fall instrument called a thermal 
gradiometer made microstructure measurements in Mc San Diego Trough, 
Florida current and in the central Pacific between San Diego and Hawaii. 
He found regions of high gradient to be about 1 meter thick and have a 


gradient around 10 times the mean gradient. Figure 8 is a record of the 


results he obtained with the gradiometer. ‘By tracing a region of high | 
“the region is at a depth greater than 400 meters. Osborn suggests that 


the microstructure is on a scale of centimeter order. 





Neshyba, Neal, and Denner [1969] made measurements of the temper- 
ature microstructure in the Arctic Ocean using an XBT lowered mechan- 
ically from Fletcher's Ice Island (T-3). The microstructure found had 
a step-like structure with layers that were extremely homogeneous as 
can be seen in Figure 9. The transistion layer was found to be less than 
20 centimeters thick. Figure 9b shows a section of the trace recorded 
at high gain (the lip at the end of a transition layer is due to recorder 

overshoot). The steps from 220 to 300 meters were found to be very 
stable. Figure 10 is a time series taken with the XBT modified to fall 
at 7 ft/sec instead of the normal 20 ft/sec, and demonstrates the stabil- 
ity. Below 350 meters, the homogeneous layers are larger, but still are 


separated by steps of temperature change. This is also apparent from 


Figure 10. It should be noted that these XBT recordings are independent 


of ship motion since the ice island is a very stable platfcrm. 
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Figure 8, Corrected data from the upper portion of the 
water column in the San Diego Trough on 18 September 1968. 
The curve on the right is the continuation of the curve on 
the left [Osborn, 1969]. 
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Figure 9. Vertical profile of temperature under T-3 
84 38'N, 128 21.6'W, 19 March 1969). Shaded section of 
profile (a) shown as observed in profile (b) to left [Neshyba, 
Neal, and Denner, 1969]. 
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Simpson [1971], using a Bathysonde CTD system, examined the density 
and thermal microstructure in the western Irish Sea. He found vertical 
layers up to 10 meters thick and in cases, less than 1 meter thick. The 
average gradient was as low as 0.002°C/meter in layers up to 10 meters 
thick and as high as 2°C/meter in the layers of sharp gradient. Figure 11 
gives three temperature profiles of the region. Horizontally the layers 


extended up to a distance of 3 kilometers, 
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Figure 11. Unretouched Bathsonde profiles from the western 
Irish Sea of temperature and conductivity showing pronounced 
microstructure [Simpson, 1971]. 


Be MICROSTRUCTURE MECHANISMS 


The exact nature of the processes producing, maintaining, or destroy- 


ing the thermal microstructure still are unclear because of the complexity 
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of the problem. Some mechanisms have been proposed to explain the xe- 
1, Insolation | 
Insolation (radiation from the sun and atmosphere) results in the 


heating of the water's surface. The heat is transported downward mainly 


by mixing caused by wave action. "The awind-wavevactionvends to.climinate 
‘the surface layer microstructure. However , Sagar [1960] suggests that 


after sunrise the insolation may be large enough to dominate over wave 
action, and so a layered structure could be transported down the thermal 
gradient. The depth that the layers reach is dependent on wave action. 
Waves are required to produce the mixing, but if they are too strong, 
then vigorous mixing will eliminate the layered structure. Sagar [1960] 
the microstructure was absent. 

2. Double Diffusion 


This mechanism results from the differences in the molecular 
diffusivity of temperature and salinity which are related by an approxi- 
mate 100:l ratio. Due to “iy differing diffusivity in regions where a 
warm saline mass of water overlies a cooler mass, the masses become 
unstable, forming a layered structure. Much evidence from ocean obser- 
vations [Cooper and Stommel, 1968, and Tait and Howe, 1968], a theor- 
etical study [Turner, 1968], and laboratory experiments [Turner and 


Stommel, 1964, Turner, 1967, and Stern, 1969] support this mechanism. 
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3. Internal Waves 
Woods [1968a] observed internal waves along an individual micro- 
structure sheet to be in long coherent trains. If these waves become 
unstable , they would break, giving up their energy to turbulent processes 
cauSing the sheets to be eliminated by mixing. Orlanski and Bryan [1969] 
have proposed that internal waves breaking could also establish a step struc~ 
ture by leaving a "blob" of water behind that would spread out into a thin 
layer. 
4, Current Shear 
Current shear can cause vertical turbulent mixing when the 
Richardson Number is iess than 0. 25 {Miles and Howard, 1963j. This 
turbulence occurs as a patch which is smeared out, leaving a scar that 
can lead to microstructure. 
Lovett [1969] and Pingre [1969] have proposed that these isolated 
turbulent patches can also be spread horizontally by turbulence and thus 
form a sheet. Dye studies illustrate this layering in the deep ocean 


[Schuert, 1969]. 


The complexity of the problem is apparent when it is realized 
that not one, but probably a combination of these mechanisms is occurring 


Simultaneously. 
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IV. INFORMATION CONCERNING THE THERMAL MICROSTRUCTURE 
iA Can BE CATHERED FROM BXISTING EXPENDABLE 


bowl EH ERMOGCRAPH (XBT) RECORDS 


env eteABILE DATA 


Introduction of the XBT system provided a method relatively 
inexpensive, yet at the same time accurate, of obtaining temperature- 
depth profiles originally developed for ASW operations, the system's 
application has been expanded to include Navy research ships and merchant 
ships ened under the Fleet Numerical Weather Central (FNWC), 
Cooperative Oceanographic Observation Pregram. The commercial ship 
routes of the Pacific and Atlantic Oceans are a principal source of ocean 
temperature structure data. XBT's taken during a cruise are forwarded 
at the end of the cruise to FNWC, Monterey for digitizing and then 
archived for further use. One copy of the digitized record along with the 
original XBT traces is sent to the National Oceanographic Data Center 
for permanent storage in their archives. As of May 1971, approximately 
60,000 XBT's had been digitized by FNWC, Monterey and some 28 ,000 of 


these were from the northern Pacific area [Dale, 1971]. 


B. EXPENDABLE BATHYTHERMOGRAPH SYSTEM 
The expendable bathythermograph system currently in use 1S manu- 
factured by the Sippican Corporation of Marion, Massachusetts. The 


system is composed of the probe (the expendable part of the system), 
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a launcher, anda recorder. It is an automatic system for measuring 


and recording sea temperature as a function of depth. The specifications 


for the system are given in Table I [Sippican Corporation, 1968]: 


Temperature range 
emesis accuracy 
Depth range 

Depth accuracy 


Deployment cycle time 


Gradient resolution 


Full scale stylus travel 


TABLES! 


28°F to 96°F 
-2°C to 35°C 


+0. 4°F 
Wea C 


0 to 1500 feet 
0 to 2500 feet 


+ 2% or 15 feet, whichever 
is greater 


990 seconds for 1500 ft drop 
63% of a step change in 
temperature in 3 feet; 

95% of a step change in 

9 feet 


1 second 


The probe is released from a cannister which pays out wire to 


compensate for the movement of the launcher platform while the probe 


itself pays out wire to compensate for its fall. The probe falls with a 


velocity of 20 ft/sec, and has been found to have an average terminal 


velocity of 19.9 + 0.7 ft/sec [Gouzie, Sanders, and Littlehale, 1966]. 


A thermistor is mounted in the nose of the probe, and as the probe 


descends through the water column, the water temperature is sensed 
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as it flows past the thermistor and out the tail of the probe. The 
response time for the thermistor is 110 milliseconds nominal and 130 
milliseconds maximum [Choate, 1970]. The temperature resolution 
of the probe is dependent on the magnitude of the temperature change 
over a certain depth interval. 

The recorder uses a pressure sensitive paper for recording the 
temperature profile. The paper moves vertically at a constant rate 
proportional to the fall rate of the probe. The chart paper displays 
311,656 £t/in or 92. 307m/in and 9.714°F/in or 5, 286°C/in, giving a 


resolution of 2 ft and 0.1°F for visual redding. 


C. DIGITIZING 

The digitizing process used by FNWC, Monterey is completely 
described by Bauer [1969]. The CALMA Model 480 digitizer isa X- Y 
recorder that records at least every 0.15 inch of movement in either 
the X or Y direction. The speed the stylus is moved down the trace 
determines the number of points recorded. The maximum resolution 
is 0,01 inches of movement in the X or Y directicn which is equivalent 
to .053°C or .097°F in temperature and 0.92 meters or 3.12 feet in 
depth, Each trace is digitized twice. 

The digitized data is further processed using a DIGITBT computer 
program. The program first puts the data in a common format.- Then 
the two digitized records are compared point by point in the X and Y 


dimensions. If the maximum differences in both directions are less 
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than 0.035 inches (0.186°C or 0. 340°F in temperature and 3.22 meters 
or 10.92 feet in depth), the two records are averaged on a point by point 
basis to produce an average trace which is used for the remainder of the 
DIGITBT program. If the deviation in either direction is greater than 
0.035 inches, then the records are rejected and must be redigitized. 
The program next selects points that will be printed out and stored 
on the XBT record tapes. These points are called INFLECTION POINTS. 
Both the first and last points of the record are automatically selected 
as inflection points. The other points of the record are screened in 
"batches". For each batch the observation is found that is most differ- 
ent from that expected by linear interpolation between the end points 
of the batch. If that difference exceeds 0.015 inches (0.079°C or 
0.146°F and 1. 38 meters or 4.68 feet) then the observation is selected 
for output. Otherwise the batch expands to include the next value and 
the inflection test is repeated until an inflection is found and until the 
entire record has been screened. Any point that can be reconstructed 
by linear interpolation is not required and therefore is removed from 
the output. The DIGITBT program converts the inflection points from 
inches to meters and Centigrade degrees for output. The depth is 
rounded to the nearest whole unit and temperature to the nearest tenth 
of adegree. Figure 12 is an example output from the DIGITBT program 
and is the information stored in FNWC, Monterey's records. The title 


line contains a ship identification code, day, month, year, time, 
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Digitized record [Dale and Stevens, 1970 ]. 


latitude, longitude , base line temperature, grid type, and maximum 


depth of probe (15 for 1500 ft XBT or 25 for 2500 ft XBT). 


Figure 13 


1S a comparison of the original XBT trace and the reconstructed trace 


from a DIGITBT record similar to the example in Figure 12. 


As could 


be expected, the digitizing process has eliminated some of the small 


scale microstructure. 


For this thesis the data used is from the FNWC, Monterey's data 


bank that has been processed using the DIGITBT program. 
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Comparison of an original XBT trace with the trace 


reconstructed from the digitized data. The digitized trace is plotted 
to a scale of 90 meters/in due to computer requirements while the 


original trace is scaled to 92. 3 meters/in. 


Eigure 13. 
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D. THERMAL MICROSTRUCTURE 

Thermal microstructure for the purpose of this thesis is defined 
as che temperature fluctuations as a function of depth that exist ona 
scale that is smaller than the scale of the gross temperature structure 
which is derived from the large scale processes such as weather, latitude, 
season, and water mass transport. 

To analyze the XBT trace for microstructure, a method was sought 
to separate the temperature perturbations (microstructure) from the 
mean of the XBT trace. A method to filter the trace and subtract the 
mean portion of the curve was desired. Several type filters were consid- 


ered, A binomial filter was selected primarily because Roden [1971] had 











If the original series is denoted by s(z), then the mean series, s*(z), is 


found by: 





where w(i) 1s a smoothing array that has the following properties: 


The smoothing array is found by evaluating the below modified binomial 


sequence, 
; m! a m! P 
w(1i) = (m=1)!1! ee | 


The parameter, m, is even and controls the width (number of weighted 
points) of the array. The parameter, p, is an arbitrary exponent which 
controls the peakedness of the array. Roden found that a statistically 
stable mean series could be obtained if m = 40, with p= 0.5. Increasing 
m produces more smoothing, but the series is physically shortened as will 


be seen later. The response of the modified filter, r(k), is given by 
x(k) = 1-)  w(ijcos (ik 4 Q) 
1 


where k = wave number ( 27) wave length) and A= filtering interval 
(distance between the 41 weighted points). Figure 14 is the solution for 
the response function for filtering intervals (42) of 1, 2, and 3 meters. 
The smaller a the sharper is the response curve at a cut-off wave 
length, however, the cut-off wave length is lowered. It should be 
remembered that with a binomial filter the amplitude response 1s not very 


predictable in the region of cut-off. 
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Figure 14, Computer solution for modified binomial filter 


with filtering intervals of 1, 2, and 3 meters. 
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A theoretical XBT trace was constructed for test purposes. An 
isothermal layer of 50 meters was followed by a thermocline simulated 
with an exponential decay function (e **) that decayed in 650 meters. 
The microstructure was constructed by superimposing a sine wave 
function (sin2 7 kz) on the exponential function. The amplitude was 
varied to represent the temperature fluctuations desired, and the wave 
length varied to represent the microstructure layer thickness scale. 
Figure 15 is a composite picture of a theoretical trace, the mean 
filtered trace, and the resulting temperature perturbation curve for a 
7 meter wave length, 0.5°C temperature variation, and a filtering 
interval of I meter. As expected, the perturbation curve is the input 
Sine wave curve with both the wave length and amplitude of the temper- 
ature variation being preserved. Figure 16 is another set of theoretical 
curves showing that for the one meter Filtering interval if the wave 
length is 34 meters and the temperature amplitude variation is again 
0.5°C the temperature perturbation amplitude has been amplified, while 
the wave length is preserved, It should be noted that the 34 meter wave 
length is on the slope part of the response curve in Figure 14, and this 
shows the instability of the filter in this region as predicted. If the 


filtering interval is increased to 3 meters then the original sine wave 


amplitude will again be passed. 
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Figure 15, 
XBT trace. 
and Temperature Perturbation curve for a 7 meter wave length, 0. 5°C 
Temperature variation, and a 1 meter filtering interval. 
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Result of modified binomial filter applied to Theorectical 
From left to right, Theoretical XBT trace, Mean curve, 
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Figure 16. Result of modified binomial filter applied to theoretical 
XBT trace showing the amplification in the temperature perturbation 
amplitude when the wave length of the perturbation is on the slope 
part of the response curve of Figure 14. From left to right., Theo- 
retical XBT trace, Mean curve, and Temperature Perturbation curve 
for a 34 meter wave length, 0.5°C Temperature variation, and a 
1 meter filtering interval. 





= DATA SELECTED FOR ANALYSIS 

A selection of 129 XBT traces has been analyzed for this thesis using 
Roden's modified binomial filter with 1, 2, and 3 meter digitizing inter- 
vals. All of the traces were from the northern Pacific Ocean. Digitized 
information was furnished by the Fleet Numerical Weather Central (FNWC), 
Monterey, and was processed using the IBM 360/67 at the Naval Post- 
graduate School's W. R. Church Computer Center. The traces were 
extracted from the FNWC data bank by selecting geographic grid squares. 
The selection of these grid squares was based on three criteria. The first 
and most important was to get a sampling from the Subtropic, Subarctic, 
Kuroshio Transport and Region of Confluence climatic regions proposed by 
Tully [1964]. Secondly it was desired to secure data for comparison with 
copies of original XBT traces obtained from J.F.T. Saur of NOAA's 
National Marine Fisheries Service at the Fishery-Oceanography Center in 
La Jolla, California. The third criteria was to request a reasonable quan- 
tity of data from FNWC. Data was desired for all months of the year in 
each grid square, but unfortunately the coverage fell short of expectations. 
Another shortcoming in the data is that the Subarctic region has not been 
sampled as often with XBT as the Subtropic region and therefore, data is 


not very plentiful from this area. 
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F. DESCRIPTION OF MICROSTRUCTURE ON XBT TRACE 

Figure 17 is a composite diagram of an actual XBT trace, the 
digitized trace, a mean curve found with Roden's filter, and the temper- 
ature perturbation curve using a 1 meter filtering interval. The scale of 
the actual XBT trace is 92.3 meters/in while the other three traces are 
scaled 90 meters/in for computer requirements. By comparing the actual 
XBT trace and the digitized reproduction it can be seen that some of the 
very fine scale temperature microstructure has been eliminated by the 
digitizing process and some larger scales appear to have been amplified 
by the process. The exactness of the digitized trace is considered to be 
sufficiently reliable to achieve results that are realistic and comparable. 

The first and last 20 meters of the mean trace are lost because the 
filter used 20 points at 1 meter intervals on both sides of a depth to 
compute the mean value for that depth. If a 2 meter filtering interval 
is used, then a total of 80 meters is lost and the 3 meter interval reduces 
the length by 120 meters. Since the XBT record is comparatively short 
to begin with, a 750 meter trace would be reduced by 16 percent when using 
the 3 meter interval, and if the trace is a 1500 ft XBT (460 meters) 26 
percent of the record will be lost. 

By closely comparing the digitized XBT trace and the temperature 
perturbation curve in Figure 17, it can be seen that any time an inflection 
point (as defined in Section IV C) is encountered on the XBT trace, a change 


in the slope results on the temperature perturbation curve. The number of 
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Figure 17. Composite diagram of the original XBT trace, the 
reconstructed Digitized trace, the Mean curve found by applying 
Roden's filter, and the Temperature Perturbation curve. The 
original trace is scaled 92.3 meters/in while the remaining curves 
are scaled 90.0 meters/in. 





slope changes can be counted and divided into the length of the trace 


_ is representative of the microstructure mean vertical scale. The same 


value can be obtained by dividing by the number of inflection points that 
was selected for the digitized record. To be sure that only the micro- 
structure is being considered, the inflection points or slope changes 

above the seasonal thermocline must not be used nor the depth interval 
from the surface to the seasonal thermocline. This must be determined 
for each trace independently. The amplitude of the temperature pertur- 
bations is defined as the magnitude of the temperature difference between 
the mean curve and the XBT trace. 

Fisure 28 demonstrates how the temperature perturbations are 
affected by cnanging the filtering interval. Examination of the three 
traces below 100 meters shows that if the filtering interval is increased, 
then the shape of the perturbation curve is changed. In general, it was 
found that the mean temperature perturbation vertical thickness scale 
was increased when the filtering interval ( 12) was increased, and the 
temperature amplitude was increased, but not at all depths. It should 
also be noted that more of the trace is influenced by the mixed layer depth 
as the interval is increased, so that for the 3 meter curve 150 meters 
are not useable in computing the mean temperature perturbation vertical 


thickness scale. 
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Figure 18. 
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Previous observations of the microstructure indicate scales below a 
few meters are important. Since increasing the filtering interval pro- 
duced an increase in the mean temperature perturbation vertical thickness 
scale, the perturbation amplitude and the variance, it was decided that 
the 1 meter filtering interval was more effective for processing the 


digitized XBT data. 


G. OCEANOGRAPHIC REGIONS OF THE NORTH PACIFIC OCEAN 
1. Description 

As previously mentioned, Tully [1964] divided the north Pacific 
Ocean into seven regions to be used for examining and forecasting oceano- 
graphic characteristics and behavior on a simiiar scale as the Weather 
Services have been doing for many years. In defining the regions, he 
considered volume transport, heating and cooling, precipitation and evap- 
oration, mixing (wind and convection), internal waves, and the effects 
these processes have on changing the properties of sea water. Figure 19 
is a diagram of the Pacific Ocean and the seven regions. The Tropic, 
Subtropic, Arctic, and Subarctic regions are the water mass formation 
areas where due to slow and/or zonal transport, the water mass has time 
to establish the regional properties. The normal time required to achieve 
these properties is one heating and cooling season. 

The Tropic region is bounded by Latitudes 10°S and 10°N, . Wind is 
the dominant mixing factor. The temperature is warm with a small season- 


al variation. The flow is generally zonal with a net westward transport. 
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Figure 19, Regions of the North Pacific Ocean [Tully, 1964]. 
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Figure 20. Drift in the North Pacific Ocean [Tully. 1964]. 
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The Subtropic is the largest of the regions. The flow is zonal 
and relatively slow, as is shown in Figure 20. The evaporation exceeds 
precipitation, causing an evaporation-driven convection mixing. The density 
structure is primarily a function of the temperature structure. From 
spring to autumn more heat is added during the day than is lost during 
cooling at night, so that in September the warmest surface water occurs. 
The region extends north to between Latitudes 40°N and 45°N, 

The Subarctic region is characterized by excess precipitation 
and as a result there is a distinctive salinity structure. Fresh water 
accumulates in the upper 100 meters with a halocline in the next 100 
meters where salinity increases about 1 percent. The density in this 
region is primarily a function of salinity. Mixing is due to the wind alone. 

The Pacific Arctic region also has an excess of precipitation 
and a salinity halecline. There is an ice cover in the winter and until mid- 
summer the heat added is used to melt the ice, followed by a rapid and 
shallow warming due to the long northern days. Convective mixing is 
dominant in this region. 

The Transport regions occur at the western boundaries of the 
ocean. An intensified meridianal flow transports water between the 
regions of zonal flow (Figure 20). The water does not remain in these 
transport regions long enough to fully acquire the regional characteristics. 

The Kuroshio Current is the western boundary for the Sub- 


tropic region, It transports a thick layer of warm saline water formed 
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in the Tropics near the Philippine Islands past the coast of Japan losing 
some of its heat along the way, but never becoming fully adjusted to the 
climate of the region. Its boundaries are easily identified due to large 
horizontal temperature gradients. 

The Oyashio Current is the western boundary for the Sub- 
arctic region. It has a smaller rate and volume transport than the 
Kuroshio Current. The flow is coast hugging and becomes strongest 
during the summer when it is swelled because of the ice melting in the 
PeeCciC, 

The Confluence Region is the area where the Kuroshio and 
Oyashio Currents meet. The resultant flowis eastward. The boundaries 
between these two flows can have temperature gradients as high as 5°C 
ina mile. Eddies from one stream can be dissolved in the other stream. 
As the Confluence moves east it looses its identity and is dissipated in 
the slow eastward trans-ocean flow (Figure 20). More than two years 
are required for the water to cross from Asia to America [Uda, 1963] 
so the surface processes can fully re-establish the regional characteristics. 

2. Microstructure Distribution 
The data available to this author necessitated that the regions 
investigated be restricted to the Subtropic, Subarctic, Kuroshio Current 
and Confluence regions. The first region and the one most thoroughly 


examined was the Subtropic for the purpose of comparing the results with 
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those of Roden [1971] who had conducted a detailed study of the vertical 
temperature and salinity structure in the central North Pacific in April 
1968 using STD data. 

Six locations were selected in the Subtropic region to see if 
the mean temperature perturbation vertical thickness scale was related 
within the region. The mean temperature perturbation vertical thickness 
scale was found for each of the 68 available traces using a 1 meter filter- 
ing interval. The average was then found for the traces from each of the 
six locations. This average mean temperature perturbation vertical 
thickness scale will be referred to as the average thickness scale. anes 
for each trace, the maximum amplitude of the temperature perturbation 
was found and then averaged for each location. Table II presents the 
results for the Subtropic region. The minimum, maximum, and average 
values for each location are listed along with the number of traces analyzed 
for each location. Roden's values are given for similar locations for 
comparison. The average thickness scale varies between 8.7 and 15.0 
meters which compares favorably with Roden's results which show a varia- 
tion from 13 to 16 meters for this region. One reason that a wide range 
of variations was found is that in analyzing the data, Roden used zero 
crossings to compute the thickness scale while in this analysis, inflection 
points from the digitized record were used because, as was pointed out in 
Section IV F, they represent layers resolved by the digitizing process. In 


addition, the 68 traces that were averaged represent all of the months of 
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the year and not just April. Another factor is that the maximum depth 
of the XBT traces was 750 meters while the STD record extended to 
1500 meters. 

The average maximum amplitude of the temperature perturba- 
tions varies between 0.30 to 0.49°C which is from the same minimum to 
three times the maximum value Roden found. This difference is probably 
largely due to the STD not recording the total magnitude of a temperature 
change since its response time is about one-third as fast as the XBT. 
The STD is also temperature compensated in the salinity sensor and large 
temperature changes in the water column cause a salinity spike that will 
show by comparison with the temperature record the whole change was 
not recorded. Another important factor affecting the STD is that it is 
not free from the vertical motion of the platform. This vertical motion 
as well as the physical size of the instrument could distort the temper- 
ature structure before it is sampled. 

Table III gives the results for the Subarctic region. Only 12 
traces were analyzed. The presentation is the same as for the Subtropic 
region with the exception of the omission of the Roden values since he did 
not have any results from this region. The average thickness scale varied 
between 17.4 and 82 meters. The location for the maximum thickness 
scales were in the Gulf of Alaska in the Alaskian gyre. The three traces 
analyzed were for March, April, and July with a mean thickness scale of 


84m, 101m, and 51m respectively. The general behavior (more 


37 





microstructure activity present during the long summer days and less 
during March, the coldest month) of the conditions is as would be expected 
since the area is unstable and with the cold winter,convection mixing would 
tend to homogenize the water column and with the summer heating the 
microstructure would reform. The average maximum amplitude of the 
temperature perturbations varied between 0.16 and ovante 

The Kuroshio Current and Region of Confluence results are 
shown in Table IV. The only comment on the data is that at 40°N 149°F 
the average thickness scale is smaller than for the rest of the region. 
This small value indicates a large concentration of microstructure activity 
and should be expected because this is the area where the Kuroshio and 
Oyashio currents meet, creating a region of intense vertical shear that is 
effective in producing microstructure. 

Table. V gives results from two other areas, One location was 
in the Luzon Strait between Luzon and Taiwan. This also demonstrates 
that there is a large amount of microstructure present and again is ina 
location that current shear could be expected due to the restriction of 
Flow. 

Also in Table V are the results of a 28 trace time series taken 
on two days in the California Current 200 nautical miles west of Monterey, 
California. Even though the extremes of the mean thickness scale are 
widely separated, 19 of the traces or 68 percent were within + 2 meters 


of the average value (13.0 meters). Table VI gives the local time, mean 
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thickness scale, and maximum amplitude of the temperature perturbations 
for each of the 28 traces. The hours are not consecutive because some of 
the traces were rejected by FNWC, Monterey when received for process- 
ing. Without weather information, it is not possible to correlate the 
microstructure mean thickness scales and the local time. 

The average thickness scale data from Tables 1), 1a 5 ane 
has been plotted in Figure 21 with the latitude of the location, and shows 
quite clearly that in the north Pacific Ocean for an increase in latitude 
the average mean temperature perturbation vertical thickness scale in- 
creases. A similar plot, Figure 22, shows that the average maximum 
amplitude of the temperature perturbations decreases with an increase in 
latitude. Tue to the scatter of data the relationship is net as clear as 
for the thickness scale. Figure 23 is a plot of the maximum amplitude 
of the temperature perturbations versus the average mean temperature 
perturbation vertical thickness scale. Although the scatter of data is 
wide , the relationship of the plot can be determined. For an increase in 
the thickness scale of the microstructure the maximum amplitude of the 
temperature perturbation will decrease. Sverdrup, Johnson, and Fleming 
[1942] have shown an increase in water instability when proceding north 
from Hawaii to Alaska. With an increase in instability the thick layers 
of turbulent convective mixing and of small temperature gradient would 
become thicker, eliminating the thin layers of high gradient. This is a 


possible explanation for the results shown in Figure 21. 
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Climatological data was not obtained so that no conclusions 


could be drawn concerning the local time and the microstructure activity. 
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TABLE V 


Other Areas 


California Current 36-37N , 125-126W 28 traces 


Mean Temperature MIN 
Perturbation Vertical MAX 
Thickness Scale (m) AVERAGE 


Maximum Amplitude MIN 
Of The Temperature MAX 
Perturbations (°C) AVERAGE 


Luzon Strait 20-24N , 120-124E 2 traces 


Mean Temperature MIN 
Perturbation Vertical MAX 
Thickness Scale (m) AVERAGE 


Maximum Amplitude MIN 
Of The Temperature MAX 
Perturbations (C) AVERAGE 
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TABLE VI 


California Current Time Series 


aes il Mean Temperature Maximum Amplitude 
Date mame Perturbation Vertical Of The Temperature 
Thickness Scale (m) Perturbations ?C) 
11 Nov 68 0200 15.0 0,24 
0400 15. 4 Og22 
0500 12.8 Oeo2 
0600 1343 0.16 
0700 16.6 0.20 
0900 13 0,20 
1000 hae, 0.18 
1100 P70 Ofvole 
1200 HOE Oe, 
1300 LL 0 0,25 
13 Nov 68 2300 789 0.20 
14 Nov 68 0100 eS O22 
0200 Ola? 0.16 
0400 193 0, 10 
0500 3.5 0.16 
0600 Se 0.18 
0800 L0;-7 Oe 22 
0900 iS. 2 0.20 
1000 Oe, 0. 24 
1100 1355 0.16 
1200 1229 0.14 
1300 1d, 0 0.20 
1400 10.5 0.27 
1500 eS O20 
1700 see Og22 
1800 14, 3 Uae 
1900 Zo 0. 26 
2100 12.4 0.24 
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Figure 21. Plot of the Average Mean Temperature Perturbation 
Vertical Thickness Scale versus Latitude showing an increase in 
the vertical thickness scale for an increase in latitude between 
20°N and 55°N. (4, points along Longitude 156-15PW; - , all other 
points), 
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Figure 22. Plot of the Average Maximum Amplitude of the 
Temperature Perturbations versus Latitude showing the Amplitude 
decreases with an increase in latitude between 20°N and 55°N. 

(A, points along Longitude 156-157°W; - , all other points). 
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Figure 23. Plot showing the Average Maximum Amplitude of the 
Temperature Perturbations decreases as the Average Mean Temper- 
ature Perturbation Vertical Thickness Scale increases between 
20°N and 55°N. (A, points along Longitude 156-157°W; - , all other 
points). 
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V. OPERATIONAL ANALYSIS 


A. OPERATIONAL REQUIREMENT 


Underwater acoustic energy transmissions are extremely vital to the 
Navy, especially in the anti-submarine warfare (ASW) field. Some of the 
military uses and some of the non-military uses of underwater sound as 
reported by the National Academy of Sciences [1970] are listed below by 
frequency bands. 

1. Less than 100 Hz 

a. Passive and active long range sonar 
b. SOFAR navigation and positico.ung 
c, Deep exploration and mapping 
2. 100 to 5000 Hz 
a. Shallow seismic exploration and mapping 
b. Warning bells 
c. Scuba diver communications 
d, Passive and active ship sonars 
e, Maines 
Saeco Lo KHz 
a. Underwater telephone transmissions 
b. Depth sounding equipment 
c, Fish finding 


d. High frequency ship sonar 
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46 iiseto 100 kHz 

a. Short range underwater telephone transmissions 

b. Depth soundings 

c, Fish finding 

d. Ship positioning 

e. Object location 

f. Low resolution scanning sonars 

g. High resolution sound velocity profiling 

h. Homing torpedoes 
5. Greater than 100 kHz 

a. Doppler sonar for navigation or docking 

b. High resolution scanning sonar 
For low frequency acoustic transmissions thermal microstructure 
' should not have a significant effect on the propagation because the acoustic 
wave length is large relative to the effective radius of a patch of thermal 
inhomogeneities. However, losses at low frequencies have been related to 
scattering by the index of refraction inhomogeneities [Urick, 1967]. Above 
a few hundred cycles the microstructure becomes increasingly important 
and affects sonar transmissions adversely, 

Sagar [1960], by taking temperature profiles while collecting acoustic 
(14.5 kHz) transmission data,was able to associate acoustic intensity 


fluctuations with the presence of microstructure while in the absence of 
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microstructure the variation was small and independent of the Mintzer 
dependence on range. This illustrates the necessity for the Navy to be 


concerned with thermal microstructure. 


EB. GrEererionaL METHODS. 

XBT data is being used by the Navy to obtain the gross temperature 
profile for use in the solution to ray path equations and for diagrams to 
predict sonar performance. This can be accomplished on board a ship 
or can be accomplished by Fleet Numerical Weather Central, Monterey 
in a more sophisticated manner, and the products supplied as a service to 
fleet units. 

On board ships the parameters most used are the sea surface temper- 
ature and layer depth. Tables can then be utilized to find the predicted 
range for active and/or passive sonar contacts. Little information con- 
cerning the temperature structure is presently being used in weapon 
selection or in determining the expected performance of the weapon 
once in the environment. 

FNWC, Monterey supples information on request to support ASW 
operations. The products include sonar range predictions and propagation 
loss predictions. The major factors which are used for the prediction 
computations are water mass type, bottom type, and depth of the water. 
Temperature information is furnished from the historical data bank with 


no real time input from the area in question. Historical data is used to 
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prevent the computer solution from being biased in case the ship XBT was 
from a small patch of water rather than being a representative sample of 
the total area. Again, no microstructure information is being used. 
C, POSSIBLE UTILIZATION OF MICROSTRUCTURE INFORMATION 

ON AN XBT TRACE, 

A ship has the capability of sampling the temperature in its area when 
desired. However, the XBT trace is not being fully utilized. By closer 
examination of the trace, the sonar operator could derive an indication 
of the quantity of microstructure present in his operating area, With this 
additional knowledge the operator would know whether to expect a variation 
in sonar reception or to expect steady reception. He would also know of 
the possibility of picking up contacts beyond the predicted range that 
might be held only intermittently if microstructure is present. He could 
request reduced range to maintain contact if the contact is gained at the 
predicted range and microstructure is present. Use could also be made 
of microstructure information in selecting the range for launching a homing 
torpedo. The variation in acoustic intensity will be larger for the frequen- 
cy range of the homing torpedo than it will be for sonar frequencies, so 
it would enhance the opportunity for success of the torpedo if the presence 
or absence of thermal microstructure is considered. 

FNWC, Monterey could aid the ASW problem by predicting the relative 
magnitude of microstructure to be expected for an area by examining the 


microstructure on digitized traces now in the data bank for the area. This 


could be made a part of the sonar range prediction now ave.:iable. 
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VI. CONCLUSION 


A. SUMMARY 

The motive behind the investigation of the thermal microstructure 
distribution in the north Pacific Ocean pursued in this thesis has been to 
determine what information can be obtained from XBT records and then 
used to give some prediction of the influence on sonar transmissions and 
receptions. Various techniques were explored to determine a measure of 
the thermal microstructure on an XBT trace ending with the selection of 
a modified binomial filter, suggested by Roden [1971], being preferred. 
Data was extracted from Fleet Numerical Weather Central, Monterey's 
XBT data bank, for four of Tully's [1964] climatic regions of the north 


Pacific Ocean. 


B. CONCLUSIONS 

Based on the examination of the data from the Subtropic, Subarctic, 
Kuroshio Current , and Confluence regions, the following conclusions are 
offered: 

1. The modified binomial filter is a technique that will produce a 
mean temperature perturbation vertical thickness scale for the thermal 
microstructure found on an XBT trace. 

2. The digitizing process used by Fleet Numerical Weather Genera 


Monterey, divides the trace into layers of linear tempera‘ure gradient 


He 


7 ca. 





so that the use of the digitized inflection points will give the same result 
as that obtained using inflection points from the compere perturbation 
curve. 

3. The binomial filter will also produce a measure of the amplitude 
of the temperature perturbations. 

4, The average mean temperature perturbation vertical thickness 
scale increases as latitude increases, while the average maximum ampli- 
tude of the temperature perturbations decreases as latitude increases 
in the north Pacific Ocean between 20N and 55N. 

5. The average maximum amplitude of the temperature perturbations 
decreases if the average mean thickness scale is increased in the north 
Pacific Ocean, 

6. The average mean temperature perturbation vertical thickness 
scale was found to be more dependent on the latitude than Tully's Climatic 
Regions. 

7. No association could be determined between the time of day 
or season and the mean thickness scale. 

C. RECOMMENDATIONS 

The XBT digitized data could be used for a variety of studies. Recom- 
mendations for further study enumerated below are specifically related to 
the development of the data processed in this thesis. 

1. One location should be selected for an extensive study to deter- 


mine if the XBT data will reveal a correlation with time of day and season. 
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Two likely areas in the Pacific Ocean for such a study are the San Diego 
and Hawalian operating areas used by the Navy. XBT data should be very 
plentiful from these areas. 

2. Astudy should be conducted to determine the sonar transmission 
and reception variation that can be expected from the influence of thermal 
microstructure that appears on XBT traces taken during the study. 

3. Fleet Numerical Weather Central, Monterey, California, should 
investigate the feasibility of including a microstructure quantitative pre- 
diction as a part of the Sonar Performance forecast for the active sonars. 

4, Fleet Numerical Weather Central, Monterey, California, should 
add to the identification portion of each digitized XBT trace the ci1mmatol- 
ogical data, such as wind speed, sea surface temperature, and cloud cover, 


for ease in evaluation of the trace. 


13 





ia 


tel 


LIST GiaRE BERENCES 


Bauer, R. A. , 1969, Fleet Numerical Weather Central's XBT 
Digitizing System, Computer Applications Report. 


Bergmann, P. G., 1946. ''Propagation of Radiation in a Medium 
with Random Inhomogeneities ,"" Physical Review , v. 70, 
p. 486-492. 


Choate, T. V., 1970. (Marketing Manager, Sippican Corporation), 
personal communication. 


Cohen, J. S., and Weinberg, H. , 1969. CONGRATS Ray Plot 
Program, USL Technical Memorandum 2070-110-69. 


Cooper, J. W. , and Stommel, H. , 1968. "Regularly Spaced Steps 
in the Main Thermocline near Bermuda," Journal of Geo- 


physical Research, v. 73, p. 5849-5854, 


Cooper, L. H. N., 1967. "Stratification in the Deep Ocean," 
Set ees. OW 55, p. 73590. 


Dale, D. H., 1971. (Fleet Numerical Weather Central, Monterey 
California) , personal communication. 


Dale, D. H., and Stevens, P. D. , 1970. Computer Processing 


of Expendable Bathythermograph Traces, Fleet Numerical 
Weather Central, Monterey, California,Technical Note 61. 


Denner, W. W., 1971. "The Layered Microstructure and Acoustic 
Propagation in the Arctic Ocean,'"' U. S. Navy Journal of 
Underwater Acoustics, v. 21, p. 45-51. 


Gouxie, M. W., Sanders, M. R., and Littlehale, A. D. , 1966. 


Laboratory and Engineering Evaluation of Three Models of 
Expendable Bathythermographs, USL Report 732. 


Liebermann, L., 1951. "The Effect of Temperature Inhomo- 
geneities in the Ocean on the Propagation of Sound," 


Journal of the Acoustical Society of America, v. 23, 
p. 56 3-570, 


76 





M2. 


mS. 


14, 


i. 


16, 


eae 


eS. 


69 


20. 


Pile. 


Date 


Lovett, J. R., 1968. "Vertical Temperature Gradient 
Variations Related to Current Shear and Turbulence ," 


Limnology and Oceanography, v. 13, p. 127-142. 


Magruder, P. M., 1970. Some Characteristics of Temperature 
Microstructure in the Ocean, M. S. Thesis, Naval Postgraduate 


School, Monterey, California. 


Miles, J. W., and Howard, L. N., 1963. "Note ona Hetero- 
geneous Shear Flow," Journal of Fluid Mechanics, v. 20, 
Denese 2o0. 


Mintzer, D., 1953a. "Wave Propagation in a Randomly 


Inhomogeneous Medium. I," Journal of the Acoustical 
Sogietymot America , vaw2d, p. 922-927. 


Mintzer, D., 1953b. "Wave Propagation in a Randomly 


Inhomogeneous Medium. II," Journal of the Acoustical 
Society of America, v. 25, p. 1107-1111. 


Mintzer, D., 1954. "Wave Propagation in a Randomly 
Inhomogeneous Medium. III," Journal of the Acoustical 


Suciety of America, v. 26, p. 186-190. 


National Academy of Sciences, 1970. Present and Future 
Uses of Underwater Sound, Washington, D. C. 


Neshyba, S., Neal, V. T. , and Denner, W. W., 1969. "The 
Significance of Temperature Stratification in the Arctic," 
paper presented at Symposium on Military Oceanography, 
6th, Seattle, Washington. 


Orlanski, I. , and Bryan, K., 1969. "Formation of the 
Thermocline Step Structure by Large-Amplitude Internal 


Gravity Waves," Journal of Geophysical Research, v. 74, 
Dee 75-6985, 


Ocbomne T. R., 1969. Oceanic Fine Structure, Ph. D. 
Dissertation, University of California at San Diego. 


Osborn, T. R., 1970. Report of the Symposium on Turbulence 
in the Ocean , Institute of Oceanography, University of 


British Columbia, Vancouver, British Columbia, Manuscript 
Report No. 24. 


TZ, 





23. 


24, 


Pingree, R. D., 1969. "Small-scale Structure of Temper- 
ature and Salinity Near Station Cavall," Deep-Sea 
Research, v. 16, p. 275-295, 


Roden, G. I., 1971. "Spectra of North Pacific Temper- 
ature and Salinity Perturbations in the Depth Domain," 


Journal of Physical Oceanography, v. 1, p. 25-33. 


a. "Sagar, F. H., 1960. "Acoustic Intensity Fluctuations and 


26. 


oa 


Zo. 


Zo), 


20. 


Sl, 


a2, 


See 


34, 


Temperature Microstructure in the Sea," Journal of 
the Acoustical Society of America, v. 32 Pe e121 


Schuert, E. A. , 1969. Turbulent Diffusion in the Inter- 
mediate Waters of the North Pacific Ocean, Naval 
Radiological Defense Laboratory Report TR-69-16. 


Sheehy, M. J., 1950. "Transmission of 24-kc Underwater 


Sound from a Deep Source ," Journal of the Acoustical 
Society of America, v. 22, p. 24-28, 


Simpson, J. H., 1971. "Density Stratification in the 
Western Irish Sea," Deep-Sea Research, v. 18, p. 309- 
Sie, 


Sippican Corporation, 1968. Instructions for Installation, 


Operation and Maintenance of Sippican Expendable Bathy- 
thermograph System, Sippican Report R-467B. 


Skudrzyk, E., 1957. "Scattering in an Inhomogeneous 


Medium ,"' Journal of the Acoustical Society of America, 
wee p. 50-c0. 


Stern, M. E., 1969. "Collective Instability of Salt Fingers ," 
journal of Fluid Mechanics, v. 35, p. 209-218. 


Stommel, H. , and Federov, K. N., 1967. "Small Scale 
Structure in Temperature and Salinity Near Timor and 
Mindanao," Tellus, v. 19, p. 306-325. 


Sverdrup, H. U., Johnson, M. W., and Fleming, R. H. , 1942. 
The Oceans, Prentice-Hall. 


Tait, R. I., and Howe, M. R., 1968. "Some Observations 
on Thermo-Haline Stratification in the Deep Ocean," 


Deep-Sea Research, v. 15, p. 275-280. 


78 








BO% 


36. 


37. 


So. 


Soe 


AQ, 


Al, 


42. 


43. 


Tully, J. P., 1964. "Oceanographic Regions and Assessment 
of Temperature Structure in the Seasonal Zone of the 


North Pacific Ocean," Journal of Fisheries Research 
Board of Canada, v. 21, p. 941-970. 


Turner, J. S., 1967. "Salt Fingers Across a Density 
Interface," Deep-Sea Research, v. 14, p. 599-611. 


Turner, J. S., 1968. "The Behavior of a Stable Salinity 
Gradient Heated from Below," Journal of Fluid Mechanics, 
Vso, p. 3ia-Z200. 


Turner, J. S., and Stommel, H., 1964. "A New Case of 
Convection in the Presence of Combined Vertical Salinity 
and Temperature Gradients ," Proceedings National 


Academy of Sciences, v. 52, p. 49-53. 


Uda, M., 1963. "Oceanography of the Subarctic Pacific 
Ocean," Journal of Fisheries Research Board of Canada, 
Veo pe 20-179., 





Urick, R. J. , 1967. Principles of Underwater Sound for 


Engineers , MNceGraw~Hill, p, 122-100. 
Witttmemen, D, C., Skudrzyk, E., and Urick, R. J. , 1957, 
"Forward Scattering of Sound in the Sea and Its 


Correlation with the Temperature Microstructure," 


Journal of the Acoustical Society of America, v. 29, 
p. 1124-1143. 


Woods, J. D. , 1968a. "“Wave-induced Shear Instability in 
the Summer Thermocline ,"" Journal of Fluid Mechanics 
we on, p. 79-300: 


Woods, J. D. , 1968b. "An Investigation of Some Physical 
Processes Associated with the Vertical Flow of Heat 
Through the Upper Ocean ," Meteorological Magazine, 
Ree? , De Co-7 2. 


79 





INITIAL DISTRIBUTION LIST 
No. of Copies 


1, The Oceanographer of the Navy 1 
The Madison Building 
732 N. Washington Street 
Alexandria, Virginia 22314 


2. Commanding Officer 1 
Fleet Numerical Weather Central 
Monterey, California 93940 


3. Defense Documentation Center 2 
Cameron Station 
Alexandria, Virginia 22314 


4, Library, Code 0212 J) 
Naval Postgraduate School 
Monterey, California 93940 


5. Professor W. W. Denner, Code 58Dw 5 
Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 


6. Professor N. E. J. Boston, Code 58Bb 1 
Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 


7. Department of Oceanography 3 
Naval Postgraduate School 
Monterey, California 93940 


Oo eeewR EF, M. Hunt.,. Jr., USN 5 
Wiss Corry (DD817) 
Fleet Post Office 
New York, New York 09501 


9. Dr. Ned A. Ostenso 1 
Office of Naval Research 
Code 480D 
Arlington, Virginia 22217 


80 





TUL 


1a 


L. 


Se 


Thine, J. E. 1. Saur 

U.S. Department of Commerce 

National Oceanic and Atmospheric Administration 
National Marine Fisheries Service 
Fishery-Oceanography Center 

P.O, (Bexr271 

La Jolla, California 92037 


Mr. Dean Dale 
Fleet Numerical Weather Central 
Monterey, California 93940 


Mr. Paul Stevens 
Fleet Numerical Weather Central 
Monterey, California 93940 


Mize J. Ree Claxk 
Fleet Numerical Weather Central 
Monterey, California 93940 


No. of Copies 


1 





Unclassitied 
Securtty Classification 


DOCUMENT CONTROL DATA-R& 


iSecurity classification of title, bodv of abstract and indexing annotation must be entered when the overall report is classified) 


2a. REPORT SECURITY CLASSIFICATION ; 
Unclassified 
REPORT TITLE 


An Examination of Microthermal Structure Statistics as Calculated from 
Expendable Bathythermograph Records 








ORIGINATING ACTIVITY (Corporate author) 


Naval Postgraduate School 
Monterey, California 93940 






4 OESCRIPTIVE NOTES (Type of report and inctusive dates) 


Master's Thesis; (September 1971) 


$ AUTHORIS! (First name, middie initial, last name) 


Frank M. Hunt, Jr. 


6 REPORT OATE 7a. TOTAL NO. OF PAGES 7b. NO. OF REFS 
| September 1971 83 43 


/ 8a. CONTRACT OR GRANT NO 8@. ORIGINATOR’S REPORT NUMBER(S) 


b. PROVECT NO. 





pee 9b. OTHER REPORT NO(S) (Any cther numbers that may be aasigned 
this report) 


10. OISTRIBUTION STATEMENT 


Approved for public release; distribution unlimited 











- SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 
Naval Postgraduate School 
Monterey, California 93940 





13. ABSTRACT 


Expendable Bathythermograph (XBT) temperature depth profile data from 17 locations 
jin the north Pacific Ocean has been analyzed to extract information concerning the thermal 
microstructure, The microstructure information is then examined to determine if the 
distribution can be associated with the climatic regions of the north Pacific Ocean pro- 
posed by Tully [1964]. 


A modified binomial filter used by Roden [1971] to analyze STD data was applied to 
find a mean temperature series. A temperature perturbation series results when the 
mean is subtracted from the original XBT series. Two resulting parameters, average 
mean temperature perturbation vertical thickness scale and average maximum amplitude 
of the temperature perturbations, can be associated with the microstructure. 


The average mean temperature perturbation vertical thickness scale was found to 
;increase with latitude from 20°N to 55°N, and as the thickness scale increased, the 
maximum amplitude of the perturbations was found to decrease. 





FORM (PAGE | 
DD ‘ocat & OS get \ 82 Unclaiscataied 
S/N 0101-807-6811 | 


Secur: , Classification 


a 


-31408 





Unclassified 


Security Classification 


KEY WORDS 


Microstructure 
Microthermal structure 
Expendable Bathythermograph 
XBT 

Digitizing 

Temperature 

Underwater Sound 

Climatic Regions 


Temperature Perturbations 


9 D i 4/3 (BACK ) Unclass: ied 


JN 0101-°807-6823 83 a 





Security Clas .fication A-39 4 4 




















7 JUL 73 2 lec Bald 
| 16 OCT 74 ec7su 
| 30 JUL 78 24/42 





Thesis 130726 


H9420 Hunt 

e.2 An examination of 
microthermal structure 
Statistics as calcula- 
ted from expendable 
bathythermograph re- 
cords. 


‘il il lI | il i 


f Hl i 


oo et timer 





